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ENGINE ANGULAR MOMENTUM ON LONGITUDINAL AND DIRECTIONAL 
STABILITY IN STEADY ROLLING MANEUVERS 1 
By Ordway B. Gates, Jr., and C. H. Wood H 1. n g 


SUMMARY 


The effect of engine momentum on the longitudinal and directional 
stability of aircraft in steady rolling maneuvers has been investigated. 
The results presented indicate that the gyroscopic moments produced on 
the aircraft by a rotating engine in rolling maneuvers can have an 
appreciable effect on the range of rolling velocities for -which longi- 
tudinal or directional instability might occur. 


INTRODUCTION 


The analysis presented in reference 1 of the effect of steady rolling 
on the longitudinal and directional stability of aircraft was made for the 
assumption of zero engine momentum; hence, the results presented were 
independent of the direction of rolling. Some of the present-day aircraft 
have exhibited different characteristics in left and right rolls which can 
be attributed to the asymmetric moments produced on the aircraft by the 
rotating engine. The purpose of this analysis is to present the aircraft 
equations which include these asymmetric engine gyroscopic moments and to 
demonstrate the effects of these terms on the divergence boundaries pre- 
sented in -reference 1 for the steady rolling case. The divergence bound- 
aries presented in this paper are for aircraft having static stability. 


SYMBOLS 


W weight, lb 

T total aerodynamic moment, lb -ft 

■'"Supersedes recently declassified NACA Research Memorandum L 55 GO 5 
by Ordway B. Gates, Jr., and C. H. Woodling, 1955- 
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Mjj rolling moment^ — lb-f*t - ; 

My pit citing moment, lb -ft 

i 

M 2 yawing moment, Ib-ft 

H angular momentum, slug-ft 2 / sec 

Iy moment of inertia about-body X-axis, slug-ft 2 

Iy moment of inertia about body Y-axis, slug-ft 2 

12 moment of inertia about body Z-axis, slug-ft 2 

Iy2 product of inertia in X, Z-plane (positive when principal 

X-axis is below body X-axis at nose), slug-ffc2 

Ix e moment of- inertia of engine about body X-axis, slug-ft 2 

p rotational velocity about body X-axis, radians/sec 

q rotational velocity about body Y-axis, radians/sec 

r rotational velocity about body Z-axis, radians/sec 

03 e engine rotational velocity, radians/sec * 

V aircraft velocity, ft/sec 

u component of V along X-body axis, ft/sec 

v component of- V along Y-body axis, ft/sec 

w component of V along Z-body axis, ft/sec 

a angle of attack of X-body axis, w/u, radians 

3 angle of sideslip, v/u, radians 

F total aerodynamic force, lb 

Fx, Fy, F 2 components of- the aerodynamic force along X-, Y-, and 
Z-axes, lb 

direction cosines relating aircraft body axes to earth 

Z-axis ► 


IH5, n^ 
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t time, sec 



P 

!q, ratio of damping to critical damping in pitch and yaw, respectively 

a^, a^, a^, a^, a Q coefficients of characteristic equation 

D differential operator, d( ) /dt 

-5_ total time derivative 

Dt 

A dot over a symbol indicates differentiation with respect to time. 

EQUATIONS OF MOTION 


The aircraft equations of motion written in vector form are: 
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-£-(5) = T 



F + W J 


( 1 ) 


where refers to total differentiation with respect to time, 2 is 

Dt _ 

the angular momentum vector, X represents the aerodynamically applied 
moments, V 7 is the aircraft velocity vector, F represents the aerody- 
naxaically applied forces, and W Is the weight of the- aircraft. A 
right-handed system of axes is chosen which originates at the center of 
gravity of the aircraft and which is fixed in the aircraft. The X-axis 
is assumed to be coincident with the X-axis of the engine and the 
X,Z-plane is considered to be the plane of symmetry of the aircraft. 
Also, the mass distribution of the engine is assumed to be symmetrical 
about the X-axis, and the engine is assumed to be rotating with constant 
speed. The rotational velocities about the X-, Y-, and Z-axes are p, 
q, and r, respectively, and the components of V in this system of 
axes are u, v, and w. 


For the previous conditions, the vector momentum is given by 
H - 4- JX^q + k(l z r - I^p) 

where X, 3", and k are unit— vectors in X-, Y-, and Z-directlons; X x , 
Iy, and 1^ are the moments of inertia; %xz 1® the product of inertia 
in the XZ-plane; is the moment of inertia of the engine about the 

X-axis; and <u e is the rotational velocity of the engine about- this 
axis, taken positive in the same sense as the rolling velocity p. 

Equations (l) become, after differentiation and resolution into 
X— , Y-, and Z- components: 

Rolling : 

r XP X TZF ~ ( X Y ~ I z) ( l r " I XZP < 1 = XI^X C 2 ^) 

Pitching: 

iyq - (l z - I X )pr + *xz(p 2 - r2 ) + = XI Ny 


(2b) 
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Yawing: 


x z *■ - - I xzp + *XZ<l r - “z ( 2c > 


X- force: 


m( u + qw - vr ) = ^ Fy + 


Viz 


(2d) 


Y-force: 


a(v - pw + ur) = Fy + Wm.^ 


(2e) 


Z-force: 


m(w + pv - uq.) = / F z + Wn^ 


(2f) 


The terms Z^, m^, and are the direction cosines between the earth 

Z-axis and the axes being used, which are fixed in the body. The equa- 
tions which relate Zj, and n^ to the airplane rotational velocities 

p, q, and r are 

Z 5 = m 5 r - n 5 q 
= n^p - Z ? r 
n^ = Z ? q - m^p 

Certain assumptions beyond that of a constant rolling velocity are 
necessary in order to linearize these equations. The term p2 - in 

the pitching equation (2b) is taken as approximately equal to p2; the term 
Ixzqr in the yawing equation (2c) is considered negligible; and it is 
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assumed tha t no change occurs In the X- component of the forward velocity. 
Further, It Is assumed that 


— » p - » a 

u u 

and, for the ass ump tion of— small out-of-trim aerody nami c forces, the 
Y- force and Z-force equations (2e) and (2f) are approximately 

P = pa - r 


& = q - pP 


Also, the aerodynamic moments are taken as 

My = Mqq. + At 


M z = N r r + NpP 


Equations (2), for these assumptions and substitutions, become in determi- 
nant form - 


"“If 


Ma 

C Y 


~ (*X - I y)p q - 1 X G (a e 


An, 

*Y 


N, 


P 


(iX - Iz)Po + 


D _ — 
Iv 


Iv ° 


= 0 


0 

-1 


-P. 


D 


D 
P r 


1 

0 


= Po°o 
= 0 


where p Q refers to a constant value of the rolling velocity p and 

Oq, to the initial value of the angle of attack. Expansion of this 

determinant with the right-hand side set equal to zero yields a charac- 
teristic equation of the form 
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ai|D 


a^I^ + SgD^ + ayJD + Bq — 0 


(3) 


where 


a 4 

a 3 


= 1 


_Nr _ % 

X Z T Y 

- b 2 + 2 £_^ + ^_ 


ag — p 0 


L Z x Y 


I Y I Z 


( X Z “ ^Pq ~ I X e a - > e 


( X X ~ ^Po + I X e cu e 


*, =_^r p 2 _^k p 2 ,¥P + Mr 

1 X Z ° X Y ° X Y X Z X Y X Z 


M N r N p 

®o - tV Po 2 - P D X s 

I Y I Z ^Z 


( X Z " ^Pq " 


M a- N B _ _ ^ 

I Y I Z ° X Y 


(l x - Iy) Pq + I X e CI ^ 


( X Z " ^o ~ 


( X X - X y)po + ^X^s 


X Y 


If the following substitutions are 


“V 


2 % 


W 




2 _ 


■“a 


I Y p c 


“Mq 

**« ' ^ 


■% 
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2^ = 


- N r 


T 


0U 


P 


O 


these coefficients "become 

aij. = 1 

a 3 = p o( 2? r* + 2 ^e) 


a 2 = Po 


1 + Wy 2 + <x> 0 2 


( ] 


■z - x X " - L X e T ) ( X X ■ X Y + x X, 
I Y I Z 


- Iy - Iy. t ) f^Y - Iy + Iv T 


0 




a i = Po 5 ^/^ + 2§ e a> e + 2 £d®e“* 2 - t - 2 £/ i V 0 e 2 ) 

2 ( I z “ J x ~ I x e T ) 2 2 

— ZS—i. + m»na n + 


Sq = P, 


N-0^^ - <y 


*> 8 


m 2 

" I Y + I X e T \ 

fl Z " I X " t X^A| 

/X X “ J Y "^XgA 

“9 

* ) 1 

V ^ / 

V 12 / 


which are essentially equivalent to the coefficients presented on page 9 
of reference 1, with the exception of_the engine momentum terms. It vilL 
be noted that the p Q factors which multiply the coefficients shown 
here do not- appear in 'reference T7 This difference is attributable to the 
fact- that the differential operator used in this paper is a time operator, 
whereas in reference 1 the operator has been made a function of rolling 
velocity. 
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ANALYSIS AND DISCUSSION 

The combinations of and eng that result in a negative value 

for the constant term ao of equation ( 3 ) for given values of £ Q , 6 ^, 

and t and hence give aperiodic instability are of primary interest. 

For an airplane with given values of natural frequency in pitch and yaw, 
these values of and define the rolling velocities for which 

this instability will exist. 

A sample of these aperiodic divergence boundaries is shown in 
figure 1 for t = 0 (cue = 0) for £ 9 !^ = 0 and 0.0031* The mass and 

aerodynamic characteristics of the airplane for which these boundaries 
were constructed are presented in table I. It should be noted that 
boundaries constructed for a constant value of do not correspond 

MnN r 

to a constant value of ; instead, every point on the curve repre- 

sents, dimensionally, a different value of this parameter. In this 
plane, the frequencies in pitch and yaw of a given airplane, for all 
values of p 0 , lie along a straight line similar to the one shown in 
the figure. The point shown for p Q = 1 radian/sec defines the fre- 
quencies of the airplane chosen for this illustration. The slope of 
this line is determined from the ratio of the square of the natural 
frequencies in. pitch and yaw. For the case shown, the frequency locus 
of this airplane passes through the divergence boundary constructed for 
= 0.0031 for p c = ± 1.8 radians/ sec and remains on the unstable 

side of the boundary up to p Q = +2. 3 radians/sec. Generally, the char- 
acteristic roots of the system in the unstable region of this plane are 
a pair of stable complex roots, one stable real root, and one -unstable 
real root. 

It would be possible to take into account the engine momentum by 
plotting boundaries for various values of I^t, but, as was mentioned 

in the discussion of M q , these boundaries would not correspond to a 

constant value of engine momentum. Also, both positive and negative 
values of I x e T would have to be considered in order to cover both the 

left and right rolling conditions. The former difficulty can be avoided 
for both these cases by plotting the boundaries as a function of the 
dimensional frequency parameters rather than in terms of and cog 2 . 

Presentation in this form necessitates the construction of a boundary 
for each rolling velocity, but this construction is relatively simple. 
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MqN r 

A sample of these baumda.ri.es is Bhown in figure 2 for = 0.044 

and Xxe<% = 0. For this ease the boundaries for left and right ro ll lin g 
are identical. The effect- of the engine momentum on these boundaries 
can be seen in figure 3 for the case of Ix e <% = 17*554 slug-ft 2 /sec. 

Boundaries are presented for both right and-left rolls, and the critical 
roll velocities are shorn on the figure for both rolling conditions. 

For positive (right) rolling the unstable range of p is between 

p Q = 2.1 radians /sec and P Q = 2.5 radians/sec and for negative- ( left ) 

rolling this range is between p Q = - 1*7 radians /sec and 

P 0 = - 2.2 radians/sec. Xn order to specify the absolute range of the 

roll rate which might be critical for a given airplane, it is necessary 
to know the critical values of p for both the left and right rolls . 

For the value of considered here, this unstable range for the 

example airplane would be defined as being 


1.7 < |P 0 | < 2-5 

Curves are plotted in figure 4 to show the effect of engine momentum on 
the values of critical p (p both negative and positive) for the par- 
ticular airplane being considered. Xt can be seen that, for Xxe^te 3 0 * 
the absolute range of critical p is between Ip f of 1.8 and 2.3 radians 
per second and increases to |p| of 1.5 and 2.7 radians per second for 
Iy e <% = 40,000 slug-ft 2 /sec. Thus, the range- of ro lling velocities for 

which a given airplane might experience instability, based on this steady 
rolling assumption, can increase appreciably with the magnitude of the 
momentum of the rotating engine; hence, the effect of engine momentum 
should be considered in the analysis. 

A point of interest with respect to the construction of the diver- 
gence boundaries in the dimensional frequency plane is that it is pos- 
sible to obtain mathematically the envelope of these boundaries by plot- 
ting the locus of_ the points of maximum curvature of the respective 
curves. The mathematical expression for the curvature of a given func- 
tion can be found in any calculus’ text book (for example, ref; 2), and 
in order to obtain the desired envelope it-ls necessary only to maximize 
this expression in the proper manner ■ From this envelope and the zero 

% -Ma - 

and asymptotic values of — and , which are rather easy to calculate, 

I z I Y 
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the divergence boundaries for any given value of roll rate for specified 
inertia characteristics can be approximated with reasonable accuracy. 

In any event, the envelopes of the branches of the curves will define, 
for the steady rolling case, the combinations of pitch and yaw frequency 
for which there will be no roll- induced instability. The equations from, 
which the envelopes can be calculated, and expressions for the previously 

Hfi -M<x 

mentioned asymptotic values of — and are: 

X Y 


(a) Equations for determination of zero and asymptotic values of 

Np % 

— - and — : 

X Z X Y 


'MqX _ Po^e ~ ( X Z ~ x x)pq 2 



MqHr 2 

XfXz ° 


X X ~ I Y r 2 
X Z ° 


+ 


X X e “eP Q 

X Z 


+ 


^“ePo 

X Y 


( X Z ~ x x)Po 2 
X Y 


M -- 

( X X " x y)po 2 + x X e “ePo 

Wk. 

L iz J 

% 


(*l\ __ 

1 

CVJ 

0 

Pi 

sH? 

1 

V Iz /M ^ x X e “ePo (Xz - X X )p 0 2 

X Y 

X Y x y 




x X e “ePo 


I 


Z 
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(b ) Equations for determination of envelope of divergence boundaries 


x z 


( T X ~ X y)Pq 2 + 


fep 2 = a ± b 

% x z ° 


% 

% 


I X e a3 e^ ) o 


( X Z “ x x)Po 2 
% 



c + b 


The equations presented for determination of the envelopes of the diver- 
gence boundaries require some further explanation. The combinations of 


X Z 


and 


-Mg 

X Y 


which define the- envelopes are 


N, 

X 


a + b 


and 




% 

% 


c + b 


The branch of the f ami ly of divergence boundaries to which each combina- 
tion applies depends on the- sign of the rolling velocity. A sample of 
the envelopes calculated for the boundaries of figures 2 and 5 is shown 
in figure 5* 
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concluding remarks 


An analysis has been presented, to examine the effects of engine 
momentum on the longitudinal, and directional stability of aircraft. 

The results indicate that the ranges of rolling velocity for 'which the 
aircraft might experience a roll- induced aperiodic divergence in steady 
rolling maneuvers can be appreciably increased by the engine momentum. 
For a particular airplane used in a sample calculation, the range of 
critical rolling velocity p was calculated to be, when engine momen- 
tum was assumed zero. 


1.8 < jp] < 2.5 


and, when an engine momentum of 17 , 55 ^ slug- ft 2 / sec was considered, 
the range was extended to 


1-7 < jp{ < 2.5 


For values of engine momentum higher than that assumed for this illus- 
trative example, the range would, of course, be further expanded. 

Calculations to show the effects of including engine momentum on 
the construction of the divergence boundaries were also presented, and 
an alternate method of construction to that presented in NACA Technical 
Note 1627 was discussed. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , June JO, 1955* 
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TABLE I.- MASS AND AERODYNAMIC CHARACTERISTICS 
OF EXAMPLE AIRPLANE 


Ixj slug-ft 2 . 10,97? 

I Y > slug-ft 2 57,100 

I z , slug-ft 2 64,975 

1x2, slug-ft 2 ■ 942 

I x slug-ft 2 / sec 17,554 



1 

sec 

JL 

% sec 2 

1 

^ Z sec 2 
V, ft/sec 


-0.421 


-5.30 


2.38 

691 
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Figure 1.- Boundaries in the afy 2 , iTq 2 plane ■which define regions of 
aperiodic divergence for example aircraft. Ix^te = °* 
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P 0 ,radians/sec 



Ug_ I 
I z ’ sec 2 


(a) Right rolls. 

Figure 3 .- Boundaries in the N p /l z , M^Iy plane for example aircraft 
which define regions of aperiodic divergence as a function of rolling 

velocity. = 17,554 -- ug ~ ftg - = 0 .044. 
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P 0 , radians /sec 

1.0 1.7 2.2 2.5 



2 3 4 5 6 

M _L 

I z ’ sec 2 


(b ) Left rolls . 
Figure 3-- Concluded- 


Rolling velocity, radians/sec 



Figure 4.- Effect of engine momentum on rolling-velocity range 

for which example aircraft is unstable. 


if 
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